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Abstract. The physical modeling of the accretion disk boundary layer, the region where the 
disk meets the surface of the accreting star, usually relies on the assumption that angular mo- 
mentum transport is opposite to the radial angular frequency gradient of the disk. The standard 
' model for turbulent shear viscosity, widely adopted in astrophysics, satisfies this assumption by 
construction. However, this behavior is not supported by numerical simulations of turbulent 
magnetohydrodynamic (MHD) accretion disks, which show that angular momentum transport 
driven by the magnetorotational instability is inefficient in this inner disk region. I will discuss 
the results of a recent study on the generation of hydromagnetic stresses and energy density in 
the boundary layer around a weakly magnetized star. Our findings suggest that although mag- 
netic energy density can be significantly amplified in this region, angular momentum transport 
is rather inefficient. This seems consistent with the results obtained in numerical simulations and 
suggests that the detailed structure of turbulent MHD boundary layers could differ appreciably 
^ , from those derived within the standard framework of turbulent shear viscosity. 
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For an accretion disk to merge smoothly to the surface of the central accreting star, its 
angular velocity must decrease inwards in the inner disk to match the sub-Keplerian 
stellar spin rate fi* (jFrank et al.l 120021 : iHartmann 2009; lArinitagjlioioh . It is easy to 



derive from the standard accretion disk model that half of the accretion energy is released 
in this so called accretion disk boundary layer (ADBL). Understandi ng the detailed 



fvq , physics within this thin layer, therefore, is of great ast rophysical interest (jPiro fe Bildsten 
2004 iBalsara et ai1 l2009; llnogamov fe Sunvaevl[2010h . 



According to standard models, the structure of ADBL depends entirely on turbulent 
angular momentum transport. On the one hand, standard m odels implicitly assume tha t 
the turbulent stress is linearly proportional to the local shear ( Shakura fc Sunvaevlll973h . 
On the other hand, it is widely accepted that the magn etorotational instability (MRI ) 
is responsible to drive turbulence in the accretion disk (Ba lbus fc Hawlevlll99lill998h . 



Given that the shear rate, q = — dln51/<ilnr, is negative in the boundary layer, the MRI 
is inactive and the standard ADBL model breaks down. 

To understand the behavior of the ADBL and derive a self-consistent model, we study 
the local evolution of magnetohydrodynamic (MHD) waves when they are stable to the 
MRI. We explain how shearing MHD waves can generate strong toroidal magnetic fields 
without leading to efficient angular momentum transport. 
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2. Local Model for MHD Accretion Disk Boundary Layers 

We employ t he standard shearing b ox approximation (jrlawlev et al.lll995h and follow 



the notation in iPessah fc Chanl (|2012f) to denote the velocity and magnetic fluctuations 
by u and b, respectively. By assuming the flow is incompressible and considering only a 
single (shearing) Fourier mode at wavenumber k = {k x , fc„, fc z ), all th e non-linear terms 



in the ideal MHD equations vanish identically ( Goodman fc Xulll994 ). 

To isolate the interesting non-MRI dynamics, we focus only on z-independent pertur- 
bation. That is, we set k z = . The non-trivial MHD can be summariz ed in the following 



second-order ordinary differential equation (see lPessah fc Chanll2012l for derivation) 



^+T{r) d ^+u% x =Q. (2.1) 

<JT Z (XT 

In the above equation, the independent variable r = k x {t)/k y = qVL^t is the dimensionless 
time, where fio is the corotating frequency of the shearing box; the dependent variable b x 
is the Fourier coefficient of b x at k. The symbol u = oja/q^o is the dimensionless form of 
the (constant) Alfven frequency oja = -Bo ■ k, where Bo is the background magnetic field. 
Finally, we use the shorthand T(r) = 2r/ (r 2 + 1). Other components of the fluctuations 
can be easily derive from b x by 

u x = — (2.2) 

U! ClT 

and the divergence-less conditions V • u = V • b = 0. 

Equation (|2.1[) does not have an analytic solution. Nevertheless, it reduces to a Bessel 
equation in the limit r — > ±oo. By keeping up to the first order terms in I /lot, we 
can derive the (mean) total stress T xy = (u x u y — b x b y ) and the (mean) energy density 
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Figure 2. The filled circles represent the value of the amplification factor E+/E- for different 
values of the shear. The dashed line shows the function 10/oj 2 = 10(<?f2o/o;A) 2 , which is in good 
agreement with the numerical results for stron g shear. This as y mpto tic behavior is insensitive 
of the initial conditions as long as C_ 7^ 0. See lPessah fc Chanl ()2012l ) for details. 



E = (u 2 + b 2 )/2, where ( • ) stands for the spatial average: 

- (S*C + SC*)sin(2uJT) 

(S*C + SC*)cos(2ujt) 
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(2.3) 
(2.4) 



The complex integration constants S and C are determined by the initial conditions; and 
the star symbol denotes complex conjugation. Note that the energy balance equation 
dE/dt — qfl()T X y is satisfied up to first order. 

We plot the numerical solutions for the energy E(t) and the stress T xy (t) in Figure Q] 
by thick and thin solid lines, respectively. The dotted and dashed lines show our analytic 
approximation with two different sets of integration constants. Note that the energy in- 
creases significantly around r = although the stress keeps oscillating around zero. The 
physical meaning is that the shearing MHD waves can extract energy from the back- 
ground without increasing the time-averaged stress. Hence, the dissipation of magnetic 
energy in the ADBL is inefficient. 

Let E± = lim T _>±oo -E(r). Despite the fact that we cannot analytically solve for the 
amplification factor E + / E- , we can still draw important conclusions from the properties 
of the numerical solutions. We plot the amplification factor as a function of 1/uj in 
Figure [2j In the limit of weak shear, the pure Alfven waves cannot lead to net energy 
gain. Nevertheless, for strong shear q uja/^o, or l/ w ^ 1, the amplification factor 
follows a power law, 
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This is remarkable because the amplification is unbounded without a linear instability. 
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3. Conclusions 

Understanding the relationship between the stress and the shear rate is essential in 
modeling the ADBL. In this proceeding, we point out that the standard turbulent ac- 
cretion disk model is inapplicable in the boundary layer because the MRI is inactive 
there. We provide the key steps to solve stable shearing MHD waves, which describes 
the physics better. Although the energy of the waves can be significantly amplified, the 
time- averaged stress remains zero. Our findings agree with global MHD simu lations of 
accretion disks performed bv lArmitagel (|2002h ; ISteinacker fc Papaloizoul (|2002| ). 

Our results suggest that the ADBL will be strongly magnetized. However, the mag- 
netic fields will not be able to transport angular momentum nor dissipate to heat up 
the plasma. One direct consequence is that the angular momentum transport from an 
accretion disk to a central accreting star, as predicted by the standard ADBL models, 
is suppressed. Some other physical mechanism is therefore needed to spin up weakly 
magnetized stars. 
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